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We present a quantum Monte Carlo study of a Heisenberg antiferromagnet on a spatially anisotropic 
square lattice, where the coupling strength in the x-direction [J^) is different from that in the y- 
direction {Jy). By varying the anisotropy a from to 1, we interpolate between the one-dimensional 
chain and the two-dimensional isotropic square lattice. Both 5 = 1/2 and S = 1 systems are 
considered separately in order to facilitate comparison. The temperature dependence of the uniform 
susceptibility and the spin-spin correlation length are computed down to very low temperatures for 
various values of a. For S — 1, the existence of a quantum critical point at af = 0.040(5) as well 
as the scaling of the spin gap is confirmed. Universal quantities predicted from the 0(3) nonlinear 
G model agree with our results at a = 0.04 without any adjustable parameters. On the other hand, 
the 5=1/2 results are consistent with a^^^^'^ — 0, as discussed by a number of previous theoretical 
studies. Experimental implications for S — 1/2 compounds such as Sr2Cu03 are also discussed. 



PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

Low-dimensional quantum magnetism is currently one 
of the most intensively studied fields in condensed mat- 
ter physics. Through the synergistic efforts of theory, 
numerical simulation, and experiment, our understand- 
ing of the static and dynamic behavior of low dimen- 
sional quantum magnets has grown tremendously. The 
study of the one-dimensional (ID) quantum Heisenberg 
antiferromagnet (QHA) has a long history dating back 
to the exact solution found by Bethdil for the S = 1/2 
nearest neighbor (NN) chain. He found the ground state 
eigenfunction for this system and showed that there is 
no long-range order at T = 0, unlike its classical coun- 
terpart. The Bethe Ansatz solution, however, is peculiar 
to the S = 1/2 chain and cannot be readily generalized. 
Recent theoretical advances in quantum magnetism have 
come primarily from applications of quantum field theory 
techniques. The low-energy, long-wavelength behavior of 
quantum antiferromagnets can be mapped onto an effec- 
tive relativistic field theory. 

In his pioneering work based on the large-spin, semi- 
classical mapping of the IDQHA to the 0(3) quantum 
nonlinear a model (QNLcrM) in (1+1) dimensions^ Hal- 
dane conjectured that all half-odd-integer spin chains 
should behave qualitatively like a S* = 1/2 chain, while 
for integer spin chains the zero-temperature spin corre- 
lations should decay exponentially with distance due to 
the presence of the so-called Haldane gap, A. La two- 
dimensions, Chakravarty, Halperin, and NelsonQ have 
mapped the long-wavelength, low-temperature behavior 
of the two-dimensional (2D) QHA onto the (24-1) dimen- 
sional QNLcrM. They obtain a phase diagram for the 
system with three regimes: quantum disordered (QD), 
quantum critical (QC), and renormalized classical (RC). 
In the QD regime, the spin correlation length ^ remains 



finite even at T = 0, as in the case of the 5* = 1 spin 
chain. In the QC regime, temperature is the only rel- 
evant energy scale, and thus ^ diverges like as 
T ^ 0: T = is the quantum critical point. In the 
RC regime, the correlation length diverges exponentially 
[~ exp(l/T)], so that the ground state has long-range 
order. 

Due to the peculiar role played by the spin quan- 
tum number in the ground state properties, the IDQHA 
has been much studied since Haldane's conjecture. The 
2DQHA has also attracted considerable attention both 
for its intrinsic interest and also because the magnetism 
in the parent compounds of the high temperature super- 
conductors is well described by the S = 1/2 2DQHA. 
Recently, attention has been given to the physics of 
magnetism in intermediate dimensions; namely, systems 
which exhibit a crossover from one to two dimensions. 
One can think of two different approaches to explore 
this problem. One method is to begin with spin ladders, 
which are obtained by coupling a small number of spin 
chains, and then to increase the number of chains in the 
spin ladder until 2D physics is obtained. Alternatively, 
one can begin with an infinite number of decoupled spin 
chains, and then increase the coupling between the spin 
chains gradually until the interchain coupling becomes 
comparable to the intrachain coupling; in other words, 
one studies a spatially anisotropic square lattice quan- 
tum Heisenberg antiferromagnet (SASLQHA). Spin lad- 
ders have been studied much recently, and have revealed 
the initially surprising result that the ground state prop- 
erties of the S = l/2|Spin ladder depend on the number of 
chains in the ladder ,□ analogous to the IDQHA. Specifi- 
cally, when even numbers of chains (or legs) are coupled 
to form a ladder (even-n), they show the same univer- 
sal behavior as integer spin chains, while odd-n S* = 1/2 
ladders behave essentially like a single S" = 1/2 chain at 
low temperatures and long wavelengths. 
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Our study focuses on the SASLQHA, in which the di- 
mensional parameter (interchain coupUng) can be varied 
continuously. The Hamiltonian of the SASLQHA is es- 
sentiaUy that of parallel spin chains forming a square 
lattice with interchain coupling (J^ > 0): 



(1) 



Here i and j label lattice sites along the y-direction and 
the x-direction, respectively. We will set Jy = J and use 
a = Jx/Jy as an anisotropy parameter in this paper. We 
use units in which % = ks = QHb = ^, and set the lattice 
constant a = 1 along both the x and y directions. By 
varying < a < 1, the Hamiltonian Eq. can interpo- 
late between a single spin chain and an isotropic square 
lattice. The 5=1/2 case has drawn special attention re- 
cently, due to its conjectured relevance to the physicsi-of 
stripe structures in high temperature superconductorsJa-Q 
One of the most interesting questions in the physics of 
the SASLQHA is the nature of the ground state. Since 
the ground state of the IDQHA (a = 0) is disordered 
while the QHA on an isotropic square lattice (a = 1) has 
a Neel ordered ground state, there should be an order- 
disorder quantum transition for a critical value oia = Uc- 
For S = 1 this seems intuitively clear. As a is increased 
from zero, the Haldane gap will decrease smoothly, van- 
ishing at a non-zero a^^^- However, it is not at all obvi- 
ous what the behavior should be for = 1/2. Although it 

is widely believed tbfyt af "^^'^ = 0, there are^ few studies 
claiming otherwise.Q'Q Sakai and TakahashiS first consid- 
ered the Hamiltonian Eq. by treating the interchain 
coupling in the small a limit via mean field theory. They 
also estimated that ac~^^'^ « and af^^ ~ 0.025 by 
calculating the susceptibility via numerical Lanzcos diag- 
onalization. Azzouz and coworkers obtained essentially 
the same results from their fields-theoretical study, albeit 
with smaller af^^ ~ 0.00186Btj The renormalization- 
group argument by Affleck and coworkersllilEHl suggests 
that the answer is not universal but depends on the mi- 
croscopic details of the model. For the nearest-neighbor 
model [Eq. (|l|)], af^^^^ = was shown by zero temper- 
ature series expansions. 

In this paper we show that our quantum Monte Carlo 
study of the S* = 1/2 SASLQHA gives results which are 
consistent with the claim that, . N ji Bel psder sets in for in- 
finitesimal a ioY S = l/2.laiaMy-y. For the = 1 
SASLQHA, we show that af-i = 0.040(5) is the quan- 
tum critical point; we find that the thermodynamic prop- 
erties at this point follow the 0(3) QNLcrM predictions 
remarkably well. We compare three quantities: the di- 
mensionless ratio Sq/Txs (defined below), the temper- 
ature dependences of the correlation length ^, and the 
uniform susceptibility Xu- These quantities agree quanti- 
tatively with the QNLcrM values without any adjustable 
parameter. 

We have carried out quantum Monte Carlo simulations 
on large lattices utilizing the loop cluster algorithm. The 



lattice size has been kept at least 10 times larger than 
the calculated correlation length. The lengths and Trot- 
ter numbers of the simulated lattices are chosen so as 
to minimize any finite-size and lattice-spacing effects. 
Spin states are updated about 2^^ times to reach equilib- 
rium and then measured 2^^ times. The same algorithm 
has been applied_£ucessfully in studying spin chains and 
spin ladders .EZlil3 We compute the temperature (T) and 
anisotropy (a) dependence of the uniform susceptibility, 
X«(a,T); the correlation length, £_{a,T); the staggered 
susceptibility, Xs(a,T); and the static structure factor at 
the antiferromagnetic wave vector Q = (tTjTt), SQ{a,T) 
for both the 5' = 1/2 and S' = 1 SASLQHA. 

The structure of this paper is as follows: Our Monte 
Carlo results for the 5' = 1/2 SASLQHA are presented 
in Sec. The uniform susceptibility and the correlation 
length data are shown and discussed here. We present 
the 5 = 1 SASLQHA Monte Carlo results in Sec. ^H. 



In Sec. IV, we determine and discuss the phase diagram 
of the 5' = 1/2 and 5 = 1 SASLQHA using QNLctM 
language, especially focusing on the quantum critical be- 
havior. In Sec. implications of this study in relation to 
experiments and other related low-dimensional quantum 
magnets are discussed. 



II. RESULTS FOR S=l/2 

In Fig. |l|, we show the uniform susceptibility per spin 
for S — 1/2, as a function of T for various a in semi-log 
scale. The results for the a = 1 square lattice are taken, 
from the recent Monte Carlo study by Kim and Troyer.E3 
The crossover from ID to 2D behavior is clear in this fig- 
ure. For small a, the uniform susceptibility follows that 
of a single chain, but begins to deviate around T/ J ~ 5a. 
For larger a, Xu lies intermediate between those for a sin- 
gle chain and the square lattice, as expected. 

In the small a limit, the Monte Carlo results can be 
compared with the results of experiments on weakly cou- 
pled spin chains. In particular, Sr2Cu03 has a struc- 
ture where S = 1/2 spin chains lie along the crystal- 
lographic fe-direction. The interchain interaction along 
the a-direction is frustrated, thus making this system 
quasi two-dimensional. In the 6c-planc, the interaction 
along the c-direction is much smaller than that in the 
^-direction, therefore, Sr2Cu03 is a very good realiza- 
tion of the 5" = 1/2 SASLQHA in the small a limit. 
Specifically, the crystallographic ^-direction corresponds 
to the y-direction in the notation of the Hamiltonian, Eq. 
(^, while the c-direction corresponds to the x-direction. 
Recently Motoyama and coworkerfil measured the mag- 
netic susceptibility of Sr2Cu03. They extracted a value 
for the intrachain exchange J = 2200(200)ii' by fitting 
the result Xo the theoretical expression proposed by Eg- 
gert et al..t3 Motoyama et al. also observed the logarith- 
mic correction term predicted in the theory. However, a 
sudden drop in the susceptibility near 2QK was observed, 
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and this was attributed to the onset of three-dimensional 
(3D) Neel order. In Fig. [|, the Sr2Cu03 experimental 
results are shown as solid circles using J — Jf, — 2200^. 
As one can see from the figure, this drop can be natu- 
rally explained by a small but non-zero interchain cou- 
pling. From a comparison with our Monte Carlo results, 
one can estimate the interchain coupling in Sr2Cu03 as 
J^ = J^K 0.002 J w 0.4meV. 

A couple of points need to be mentioned in estimating 
the interchain coupling in Sr2Cu03. First, since the anal- 
ysis of the susceptibility data involves the subtraction of 
a Curie term due to impurities, subtle effects may be diffi- 
cult to interpret unambiguously. However, recent nuclear 
magnetic resonance Knight shift measurements, in which 
the Knight shift is proportional to the uniform suscep- 
tibility, sIlow complete agreement with_.the susceptibil- 
ity data.LJ Second, Schulza and WangE^I have predicted 
the staggered magnetization of this system as a func- 
tion of the interchain coupling, "t-oi-^ 0.72( Jj^/ J)^/^ from 
their interchain mean- field theory.Ea With a = 0.002, we 
obtain mo = CLp32 « 0.064/^^, which agrees with the 
experimentallyE3 determined value of 0.06(1)/xb, giving 
credence to our estimate of the interchain coupling in 
SrsCuOg. 

The spin-spin correlation length is deduced from fits of 
the calculated instantaneous spin-spin correlation func- 
tion to the asymptotic Ornstein-Zernikc (OZ) form. Only 
large distance numerical data arc included in the fits to 
ensure that the asymptotic behavior is probed. As dis- 
cussed in our-api'ious Monte Carlo studies of quantum 
spin systems jtZnia the ID OZ form works best at high 
temperatures, while the 2D OZ form works better at low 
temperatures. However, for the S=l SASLQHA we ob- 
serve a crossover to the 3D OZ form at low tempera- 
tures. This crossover of the correlation function will be 
discussed in Sec. 0. 

The correlation lengths so-obtained for S = 1/2 are 
plotted in Fig. ^jon a logarithmic scale as a function of in- 
verse temperature. Therefore, linear behavior in this plot 
corresponds to an exponential dependence on T"^, and 
the slope corresponds to the spin stiffness. At low enough 
temperatures, all data show linear behavior except for 
a = 0, signaling that a = is indeed a critical point. 
The solid lines are the results of fits to the crossover form 
suggested by Castro Neto and Honejj which interpolates 
between the IowtX expression calculated by Hasenfratz 
and Niedermayeii23 for the (2+1) dimensional QNLtrM 
and 1^ ~ at high temperature: 



e = Aexp(27rp,(a)/r)/ 



1 + :::: 



1 T 



(2) 



2 27rp,(a) 

Two adjustable parameters are used in the fitting: A 
and Ps- The anisotropy dependence of Ps(a) is used to 
provide a crossover temperature scale between quantum 
critical and renormalized classical behavior in Sec. IV. 



We find that the Ps{ct) values extracted from fitting 
and agree within the errors, as evidenced by the iden- 
tical slopes in Fig. 0. 



III. RESULTS FOR S=l 

The uniform susceptibility data for S — 1 are shown in 
Fig. H. The results for the square lattice (a = 1) a]fe.taken 
from recent QMC work by Harada and coworkers .CJ As in 
the 5=1/2 case, for small a any deviation from the sin- 
gle chain data is only observed at very low temperatures. 
As is evident in Fig. |^, one observes distinctively different 
behavior for a — 0.02 compared with that for a = 0.05. 
Specifically, Xuict = 0.02) drops to zero at low tempera- 
tures, signaling the opening of a spin gap and concomi- 
tantly a spin liquid ground state. Indeed, one can fit the 
data with the asymptotic form XuiT) ~ exp(— A/T) at 
low temperatures. The result of the fit is shown as the 
sohd hue for a = 0.02 in Fig. |. The spin gap values 
from these fits are as follows: A(q; = 0.01)/J — 0.32(1), 
A{a = 0.02)/J = 0.25(1), and A(a = 0.03)/J = 0.17(1). 
Note that the a = 0.01 and a = 0.03 data are not shown 
in Fig. ^ for graphical purposes. 

In Fig. ^ the correlation length data in the y-direction 
are plotted for S = 1. Since most correlation length 
data along the x-direction are smaller than one lattice 
constant, they are not plotted. Note the linear scale 
in this figure, unlike the logarithmic scale in Fig. |^ for 
S — 1/2. This clearly shows the 1/T dependence of 
the correlation length for a = 0.04, which thus identifies 
a = 0.04 as a quantum critical point. For smaller val- 
ues of a the correlation length saturates at low temper- 
atures: Co(a = 0.01) = 7.0(2), ^o{a = 0.02) = 10.8(5), 
^o{a = 0.03) = 19.0(5). For a > ac, the correlation 
length diverges exponentially in 1/T. Therefore, one can 
fit the data to Eq. (H) to obtain the spin stiffness Ps{a). 

Some precautions are necessary in extracting the cor- 
relation length from the correlation function via fits to 
the OZ form. The OZ form for the correlation function 
in general is 



C{r) 



A = 



(d-1) 



(3) 



where d is the dimensionality. In quantum systems at 
low temperatures, d should be replaced by the effective 
dimensionality {d + 1). Since the low-temperature and 
long-wavelength behavior of the S = 1 SASLQHA with 
non-zero a is mapped onto the (2+1) dimensional 0(3) 
QNLctM , the functional form of the correlation function 
at low temperatures and long distances should be the 
OZ form with A = 1. Similarly, the ID S = 1 spin chain 
(SASLQHA with a =r-Q) has a correlation function of 
the A = 0.5 OZ form.ta We indeed observe a crossover 
in the functional form of C(r) from A = 0.5toA = lby 
increasing a from zero to a small but non-zero value in 
the = 1 SASLQHA. In Fig. ||, the correlation function 
at the low temperature T/ J = 0.04 is plotted to illus- 
trate the difference in A. In order to show the subtle A 
dependence, 1/C(r) is multiplied by the exponential fac- 
tor, and only the e^^^^ /£,/C{r) ^ part is shown as 
a function of r/^. The A = 1 behavior is apparent for 
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a = 0.025, while A = 0.5 describes the a = data bet- 
ter. For comparison purposes, C(r) for an 5 = 1/2 chain 
(a — 0) is also plotted; in that case the A = behavior 
is quite clear. 

One should note that this crossover is observed only in 
the QD phase, since the finite correlation lengths in this 
phase allow the QMC technique to probe the low tem- 
perature behavior. On the other hand, in phases with 
a diverging correlation length, it is difficult to study the 
ground state properties with the QMC method. Con- 
sidering that the SASLQHA behaves classically at high 
temperatures, it is possible that the 5 = 1/2 data in Fig. 
H (T/J — 0.04) have not reached low enough tempera- 
tures to show the true ground state behavior of A = 0.5. 
In deducing the correlation length data shown in Fig. ^, 
the A = 1 form was used for a > 0.02 at low tempera- 
tures, while A = 0.5 was used otherwise. 

IV. QUANTUM CRITICAL POINT 

The phase diagram of the S —\ SASLQHA is sketched 
in Fig. H, where we use the energy gap A and spin stiffness 
27rps, obtained from the fits discussed in previous sections 
as the crossover energy scales. We use the terminology 
of the QNLctM to identify the different phases in 

the figure.cl The inverse anisotropy ratio 1/a corresponds 
to the couphng constant g. The 5=1/2 phase diagram 
is almost identical to that of Ref. and hence is not 
shown here. 

At the quantum phase transition in (2+1) dimensions, 
Ps and A obey the scaling laws 

Ps - {9c " gr (4) 

A~(.9-5c)^ (5) 

with a critical exponent v ~ 0.69 for the 2DQHAi] bmce 
plays the role of the coupling constant g, we plot 
Ao(q; — as a shaded line in Fig. ^. Note that Aq 
is fixed from the 5 = 1 chain value, so that there is 
no adjustable parameter. The power law scaling fits the 
gap data quite well, confirming that the low-temperature 
long- wavelength behavior of the 5 = 1 SASLQHA is con- 
sistent with that of the 0(3) QNLctM. 

To illustrate that ac — 0.04 is indeed the quantum 
critical point of the 5 = 1 SASLQHA, the dimensionless 
ratio Sq/Txs is plotted in Fig. 0(a). According to the 
QC scafing prediction for the 0(3) QNLcrM,y this ratio 
should exhibit universal behavior in the QC regime with 
the specific value Sq/Txs = 1.10(2). As shown in Fig. 
g(a), the a = 0.04 data for Sq/Txs are constant around 
1.1 at low temperatures, in quantitative agreement with 
the QC scaling value. This ratio has been sucessfuUy 
used in identifying quantum critical behavior in difSerent 
systems, for example, weakly coupled spin ladders. E3 We 
also show the correlation length multiplied by temper- 
ature in Fig. |^(b), to emphasize the 1/T dependence of 



^j, at the quantum critical point, ac = CLQ4, at low tem- 
peratures. Using the QNLaM predictionEl = c/1.04, 
we obtain the spin wave velocity of Cy/ J — 2.6, which is 
slightly larger than the ID value (2.49), but smaller than 
the 2D value (3.067). 

Another quantitative prediction from the 0(3) 
QNLcrM is the uniform susceptibilitv_In the QC regime 
the uniform susceptibility is given asEj 

Xu{T) = f^i(oo)4T, (6) 

where i7i(cx3) = 0.26(1) is a universal constant.0 Since 
the spin wave velocity is anisotropic in this case, one 
should presumably use c^Cy instead of (? in Eq. (||). 
The inset of Fig. |^ clearly shows the linear dependence 
of Xu on temperature. The fitted value of the slope is 
0.34(1). Since we do not know the value of c^,, we use 
the fitted value of the slope and Eq. (^) to estimate 
Cx- If we take the spin wave velocity value Cy deter- 
mined above, the spin wave velocity in the x-direction is 
CxjJ ~ 0.3. This gives the ratio of the spin wave veloc- 
ity Cx/cy « 0.12, which is veiSL, close to the value given, 
by the series expansion studyiiil and spin-wave theory:Ll 
Cx/cy « 0.14. One should note that the anisotropy in the 
spin-wave velocity is enhanced from the mean field value 
{cx/cy = \f(x = 0.2) due to quantum fluctuations. 

V. DISCUSSION 

The most surprising result for the 5 = 1 SASLQHA 
is the smallness of the value af^"'^ ~ 0.04, compared 
with the Haldane gap value of the 5 = 1 spin chain 
[Ah /J ~ 0.41). One would naively think that the in- 
terchain coupling should be comparable to the Haldane 
gap to overcome the large energy gap and drive the sys- 
tem to long-range order. Indeed, such heuristic reasoning 
works well in the case of coupled even-legged spin ladders, 
where the value of the spin gap and the critical interlad- 
der coupling values are ~ 0.5J and ~ 0.3 J, respectively, 
for an array of two-leg ladders. Corresponding values for 
an array of four-leg ladders are ~ O.IJ and ^ 0.07J. 

Recently there has been a number of studies on the na- 
ture of the ground state of 5 = 1/2 spin ladders. Specif- 
ically, much theoretical interest has focused on the ques- 
tion of whether or not the ground state of the 5=1/2 
antiferromagnetic two-leg ladder is the same as that of 
the 5 = 1 chain. Since the 5=1 chain is equivalent to 
the 5 = 1/2 two- leg ladder with an infinite ferromagnetic 
coupling, the above question can be rephrased as whether 
the 5=1/2 ladder goes through a phase transition when 
the inter-chain coupling is varied from a positive (anti- 
ferroniagnetic) value to a negative (ferromagnetic) value. 
WhiteEfl has shown that the antiferromagnetic ladder can 
be transformed continously to an 5 = 1 chain by switch- 
ing on an irrelevant next nearest neighbor coupling, thus 
claiming the equivalence of the two. Kim and coworkera23 
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recently have studied a two-leg ladder with various model 
parameters by a bosonization method and have provided 
evidence that the antiferromagnetic ladder and S = 1 
chain belong to different universality classes, each hav- 
ing distinct topological string order. Our results appear 
to show the distinctively different nature of the ground 
state of the S = 1 Haldane chain and that of the an- 
tiferromagnetic ladder. However, more theoretical and 
numerical studies are necessary to resolve this issue un- 
ambiguously. 

Although both the SASLQHA and spin ladders show 
crossovers between ID and 2D behavior, the intrinsic dif- 
ference between the two approaches should be noted. At 
high temperatures, the SASLQHA behaves as an array 
of decoupled chains. The 2D behavior is only observed 
at low temperatures, as shown in Fig. ^ and Fig. |^. How- 
ever, the spin ladders at high temperatures show essen- 
tially 2D physics, since the spin-spin correlation length 
is shorter than the width of the spin ladder. Therefore, 
the crossover for the spin ladders (2D ^ ID as T ^ 0) 
is complementary to that of the SASLQHA. 

In Sec. g we estimated the exchange interaction be- 
tween chains in Sr2Cu03 as ~ 0.4 meV. Greven and Bir- 
geneau have |-Gpnsidered a similar exchange interaction 
for SrCu203.E3 Namely, they have noted that the the ef- 
fective coupling between copper spins in different ladder 
planes is mediated by Sr^+ ions and they have argued 
that this interladder exchange should be about 10 meV. 
This value is also used to describe the. behavior of the 
three leg ladder compound Sr2Cu305.EZl This represents 
quite a large discrepancy between values estimated for 
seemingly similar superexchange interactions. However, 
one should observe that the geometry of the inter-plane 
coupling in the spin ladder compound SrCu203 and that 
of the inter-chain coupling in the spin chain compound 
Sr2Cu03 are quite different. Specifically, in Sr2Cu03 , 
all Cu^"'"-Sr^+-Cu^"'" bond angles are close to 90°, while 
there exist hnear Cu^"''-Sr^+-Cu^+ bonds in SrCu203 
and Sr2Cu3 05. If one assumes that the coupling between 
these Cu^^ ions is primarily due to the superexchange 
mediated by Sr^+ ions, one can infer that the effective 
inter-plane coupling in SrCu203 is mostly due to the 180° 
Cu^"'"-Sr^+-Cu^+ superexchange interaction. This leads 
us to speculate that the small value (~ 0.4 meV) of the 
interchain coupling in Sr2Cu03 is due to the absence of 
linear Cu^'*'-Sr^+-Cu^"'" superexchange paths. 

This realization is important in understanding the spin 
fluctuations in superconducting YBa2Cu306+5, since the 
bilayer coupling path in YBa2Cu306+5 is similar to the 
inter-plan&,coupling path in the spin ladder compound 
SrCu203.t3 Although the detailed nature of the hop- 
ping integral between bilayers is not known, Andersen 
et al. have shown that the hopping mediated by the yt- 
trium ioMontributes significantly to the bilayer exchange 
couplingE2l. Obviously, more theoretical calculations of 
the bilayer exchange interaction are needed. 

In summary, we have studied the S* = 1/2 and 5* — 1 
SASLQHA with the quantum Monte Carlo method. By 



varying the anisotropy a from to 1, we go continuously 
from the one-dimensional chain to the two-dimensional 
isotropic square lattice. The temperature dependence 
of the uniform susceptibility and the spin-spin correla- 
tion length are presented for various values of a. For 
5 = 1, we show that there exists a quantum critical 
point at af^^ — 0.040(5), in agreement with other an- 
alytic predictions. The power law behavior of the spin 
gap is also confirmed. In addition, universal quantities 
predicted from the QNLcrM, such as Sq/Txs « 1.1 and 
17i(oo) « 0.26, agree with our results at a = 0.04, with- 
out any adjustable parameter. The S* = 1/2 results are 

S—l/2 

consistent with ac = 0, as discussed by a number 
of previous theoretical studies. We also estimate the in- 
terchain coupling in the 5 = 1/2 compound Sr2Cu03 by 
comparing the measured uniform susceptibility data with 
our QMC results. 
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FIG. 1. The uniform susceptibility per spin, Xu for 
S = 1/2, is shown as a function of T for various a in 
semi-logarithmic scale, a = 1 data are taken from Ref. 21. 
The Sr2Cu03 experimental results are also shown; they are 
scaled to fit the Monte Carlo results and J = 22007^ is used to 
scale temperature. The solid line is a plot of the field theory 
result from Ref. 23. 
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FIG. 2. The correlation length for 5 = 1/2 as a function 
of J/T for various values of a in a logarithmic scale. Filled 
symbols denote the correlation length in the x-direction {^x), 
while open symbols are data. For small a {a < 0.05), 
the ^x's are smaller than one lattice constant and are not 
shown in the figure. Solid lines are fits to Eq. showing 
the exponential dependence of ^ on T~^. 
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FIG. 4. The correlation length along the y-direction, 
for S = 1 as a function of J/T for various values of a. Note 
that this plot is on a linear scale, so that the straight line 
through the a = 0.04 data signifies the 1/T dependence of 
^. Filled symbols show exponentially diverging correlation 
lengths, while open symbols show the saturation of the cor- 
relation length at low temperature for a QD ground state. 



FIG. 3. The uniform susceptibility per spin, Xu for S = 1, 
shown as a function of T for various values of a on a 
semi-logarithmic scale. The a = 1 data are taken from Ref. 
28. The solid line is the result of a fit to the asymptotic form 
Xu ~ exp(— A/r). Inset: The a = 0.04 data on a linear scale. 
The solid line is Jxu = O.SAT/J. 
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FIG. 5. Inverse correlation function multiplied by the ex- 
ponential factor e~'^^^ /^/C{r) plotted as a function of dis- 
tance r/^. The solid lines are the results of fits to the form 
~ with A fixed at the given values. All data are taken at 
T/J = 0.04. 
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FIG. 7. (a) Dimensionless ratio Sq /Txs as a function of 
inverse temperature for various values of a. The solid line is 
the prediction from the 0(3) QNLcrM (b) Correlation length 
multiplied by temperature as a function of inverse tempera- 
ture. The solid line corresponds to the value cm ~ 2.49. 



FIG. 6. Phase diagram of the S = 1 SASLQHA. The open 
squares are 2ixpsl J obtained by fitting the correlation length 
to the asymptotic form Eq. ^ , and the open triangles are the 
energy gap A/ J obtained by fitting the uniform susceptibility 
to the form ~ exp( — A/T). The shaded lines are the power 
laws, Eq. (§) and Eq. (|). 
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